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^Lj Abstract: We present results for the 0{as) virtual corrections to qg — )■ Wbbq' obtained 

^ ■ with a new automatized approach to the evaluation of one-loop amplitudes in terms of 



Feynman diagrams. Together with the 0{as) corrections to qq' — ?> Wbbg, which can be 
obtained from our results by crossing symmetry, this represents the bulk of the next- 
to-leading order virtual QCD corrections to Wbb + j and Wb + j hadronic production, 
calculated in a fixed-flavor scheme with four light flavors. Furthermore, these corrections 
represent a well defined and independent subset of the 1-loop amplitudes needed for the 
NNLO calculation of Wbb. Our approach was tested against several existing results for NLO 
amplitudes including selected 0{as) one-loop corrections to W + 3j hadronic production. 
We discuss the efficiency of our method both with respect to evaluation time and numerical 
stability. 
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1 Introduction 

One of the most important background processes to single-top production as well as 
searches for a light Higgs boson is the associated production of a W boson with a pair 
of massive bottom quarks, contributing to both the W + 6-jet and W + 2 6-jet signatures. 
The precise theoretical knowledge of these processes provides moreover an excellent probe 
of our current understanding of QCD when confronted with measurements in various kine- 
matic regimes at high-energy hadron colliders. 

At the Tevatron pp collider at Fermilab, the cross sections for W boson + h jets 
production have been measured in different forms by both the CDF [1] and DO [2] col- 
laborations and the Large Hadron pp Collider (LHC) at CERN is expected to provide 
additional experimental insight with increased precision. A first measurement from the 
ATLAS collaboration has recently been published [3] and new measurements with better 
statistics are expected by the end of the year. 

On the theory side, next-to-leading order (NLO) QCD corrections to W production 
with up to two jets containing at most one heavy h jet are known [4] and cross sections for 
W+2 6-jets were determined both in the massless 6-quark approximation [5-9] and including 
6-quark mass effects [10-13] at the same level of precision. From existing NLO QCD 
calculations, the theoretical prediction for the production of W+2 jets with at least one b jet 
has been provided in Ref. [14] and compared to the CDF [1] and ATLAS [3] measurements 
in Refs. [15] and [16] respectively. Recently, the NLO calculation of Wbb has been interfaced 
with shower Monte Carlo generators using both POWHEG [17] and MCONLO [18]. The study 
of the NLO corrections to Wbb shows, however, that the NLO theoretical prediction for 
Wbb production is plagued by large renormalization and factorization scale uncertainties in 
particular at the LHC [12]. In fact, at this order of perturbative QCD, a new qg initiated 
channel with an additional parton in the final state {qg Wbbq') opens up and, being a tree 
level process, introduces a strong scale dependence. This effect is particularly pronounced 
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at the LHC, where the NLO qg channel competes with the qq' channel due to the substantial 
initial state gluon density. Only a complete NNLO calculation of pp{pp) — > Wbb can 
be expected to reduce this spurious scale dependence and give a theoretical prediction 
consistent at this order of QCD. However, this requires the evaluation of two-loop virtual 
corrections to a massive 2 — )• 3 process as well as one-loop corrections to 2 — ?> 4 massive 
processes, and single and double particle emissions through real corrections: a truly difficult 
task. In this paper we focus on and present results for one of the many contributions to 
pp{pp) — 7> Wbb at NNLO: the 0(a<j) virtual corrections to the qg Wbbq' channel, keeping 
the full bottom-quark mass dependence. Results for the qg Wbbq' channel are identical 
at the partonic level and in the following it will be understood that qg — > Wbbq' refers to 
both channels. 

Our choice is motivated by the following considerations. First of all, the 0{as) virtual 
corrections to the qg Wbbq' channel are a well-defined independent piece of the overall 
NNLO calculation of Wbb hadroproduction. It can be directly translated into the analogous 
0{as) virtual corrections to qq' — )• Wbbg by partial crossing of initial and final states. 
Once interfered with the corresponding tree level amplitudes, they provide a self-standing 
and well-defined part of the one-loop contributions to the full NNLO Wbb cross section, 
namely the one-loop virtual contributions from 2 — )• 4 processes. The remaining one-loop 
corrections come from the interference of the one-loop amplitude for the 2 — >■ 3 process 
{qq' Wbb) with itself, and is not part of this study. Furthermore, when complemented 
with the corresponding real corrections to qg — ?> Wbbq' and qq' — Wbbg, our calculation 
completely determine the NLO cross sections for both Wbb+j and, within a fully consistent 
four-flavor-number scheme, Wb + j production, i.e. for the production of aW boson with 
one or two b jets plus a light jet, where the difference between the two processes is just 
the number of b jets tagged in the final state (the parton level processes being the same in 
the four-flavor- number scheme). Since NLO real-emission contributions nowadays can be 
determined in a mostly automatized fashion with the help of existing packages, for instance 
SHERPA [19], as well as NLO Monte Carlo frameworks as POWHEG [20, 21] and MCSNLO [22], 
the virtual one-loop corrections that we present in this paper constitute the only missing 
piece for the NLO QCD cross section prediction of the above processes and are therefore 
highly desirable. Both in this and the previous case, our calculation should contribute to 
reduce the theoretical uncertainty from the unphysical scale dependence that plague the 
prediction of + 6-jets cross sections. 

The computation of QCD one-loop corrections to 2 — t- n processes with n > 4 is, even 
in the massless limit, a challenging task and retaining the mass dependence on internal 
and external particles increases the complexity even further. Therefore, only a few full 
NLO cross section computations of 2 ^ 4 and 2—7-5 processes have been carried out to 
date, among them the productions of W + 4j [23], Z + 4j [24], 4j [25], Z/j + 3j [26], 
W + 3j [27-31], W + 2'y + j [32], tijj [34, 35], ttbb [36-39], bbbb [40], W+W^jj [41, 42] 
and W~^W^bb [43, 44]. A particularly difficult component of this kind of high-multiplicity 
processes is the calculation of virtual one-loop QCD corrections. The prospect of increas- 
ingly accurate measurements at the LHC triggered a lot of interest in the improvement 
and automatization of NLO cross section predictions. Primarily two different strategies 
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have been developed for the evaluation of one-loop corrections: the traditional Feynman- 
diagram-based approach as well as unitarity techniques [45, 46]. Powerful packages like 
BlackHat [27, 47], CutTools [48], Helac-nlo [49], Rocket [29, 50] and MadLoop [51] exist 
that provide automatization and efficient numerical implementations of unitarity methods 
and that have been successfully applied to the calculation of cutting-edge one-loop pro- 
cesses. Recently, the automatized package GoSam [52] has been developed and applied to 
the automatized computation of a wide range of NLO cross sections. Moreover, several fast 
and efficient private codes exist that follow the traditional approach of Feynman diagrams 
and tensor-integral reduction. 

In this work, we develop and describe a new automatized approach to one-loop cal- 
culations based on Feynman diagrams. We test our techniques against several 2 — > 3 and 
2 — 7> 4 processes for which results are available. For instance, we are able to reproduce the 
ud — > Wddg results for W + 3j production at NLO [26, 28]. Finally, we apply them to 
the novel computation of one-loop corrections to the cross section of qg Wbbq' with a 
massive b quark. 

The paper is organized as follows: In section 2, we present the anatomy of the parton- 
level processes and discuss the general strategy to generate and simplify amplitude-specific 
expressions. Our numerically stable approach to the evaluation of one-loop tensor integrals 
is explained in section 3. Section 4 is devoted to the analysis of the achieved accuracy and 
computation times. Finally, we conclude by presenting a numerical result for dg — > Wbbu at 
NLO for a single phase-space point in section 5. Section 6 contains some brief conclusions. 

2 General strategy 

At leading order in the strong coupling, the qg — t- Wbbq' process, with the choice q = d 
and q' = u which we consider in the following, consists of 12 tree level diagrams. Examples 
of these diagrams are depicted in fig. 1 and one-loop QCD corrections are obtained by 
adding virtual gluons and fermions, yielding 308 Feynman diagrams. Ultra-violet (UV) 
and infrared (IR) divergences are regularized with dimensional regularization in d = 4 — 2e 
dimensions and we keep the full bottom-quark mass dependence while lighter quarks are 
treated as massless. We enforce transversality of external bosons through pw ■ ew = 
and Pg ■ = 0, with and p'^/g being the polarization vectors and momenta of the 
W boson and gluon, respectively. While this choice is obvious for gluons, it is justified 
for the W boson only for weak couplings to massless fermions, which is the case in the 
amplitude at hand. 

In the limit of a vanishing bottom-quark mass, the dg — )■ Wbbu transition also con- 
tributes to W + 3j production at NLO. To verify the correctness of our approach we choose 
the W + 3j subprocess ud — > Wddg, parts of which can be obtained by crossing of the 
dg —7- Wbbu diagrams and setting the bottom-quark mass to zero. Because d quarks ap- 
pear both in the initial and final state, however, the number of diagrams is doubled both 
at LO and at NLO due to contributions like the ones depicted in fig. 2 and their one-loop 
corrections. We checked our results against Ref. [28] and found agreement. 
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Figure 1. Example of LO diagrams for dg — ?> Wbbu. 
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Figure 2. Examples of additional LO diagrams required for ud Wddg which can not be obtained 
through crossing of the massless limit of dg Wbbu. 



In the traditional Feynman diagram based approach to the evaluation of one-loop 
corrections, a given NLO amplitude TW^^^ is commonly decomposed as 



(1) 



(2.1) 



with color structures Cj and polarization/spin information . Ij denotes tensor one-loop 
integrals after decomposition into tensor structures of external momenta and contraction of 
all Lorentz indices. The sum over j in eqn. (2.1) for a given term Ci runs over all one-loop 
sub-diagrams with the same color structure. 

The color- and spin-summed and/or -averaged squared amplitude is given by 



\ colors pol / \ n / 



with 



ij colors pol 

where M^^^ is the leading-order amplitude 



(2.2) 



(2.3) 



(2.4) 



decomposed in color space on the same basis of color structures Ci. 
After organizing the NLO amplitude by color factor, standard SU{3) relations are ap- 
plied to simplify the color structures and the resulting set of color coefficients is extracted. 
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Summing/averaging over final/initial color indices, after contraction with the LO color 
components, yielding the CiC* term in eqn. (2.3), is performed at this point. In the next 
step, tensor integrals are decomposed into Lorentz invariant tensor coefficients and a stan- 
dard ordering of Dirac and spinor structures is achieved with the help of anti-commutation 
relations of Dirac matrices and the application of the equations of motion. The amplitude 
is subsequently expanded in {d — 4) and UV/IR divergences are separated such that four- 
dimensional identities can be safely used without introducing the need for rational terms 
of either IR or UV origin [36]. Moreover, this approach also avoids ambiguities in the defi- 
nition of the 75 matrix, which we treat in naive dimensional regularization. The complete 
polarization information of the amplitude is contained in Dirac chains and polarization 
vectors of external bosons, commonly called standard matrix elements (SME), A4k- At 
this stage, the number of SME is of the order of several thousand for both 2 — )■ 4 processes. 
Reducing the set of SME to linear combinations in a smaller basis {TWfc} is crucial since 
the size of final expressions, and therefore the computational complexity, scales with the 
number of SME. Algebraic relations based on four-dimensional identities tailored for the 
specific process and SME at hand have been described in [36-38, 53, 54] and successfully 
applied in several calculations. However, a reduction to a sufficiently small basis is not 
straightforward, requires careful inspection of the individual contributions on a case by 
case basis, and is also dependent on the order of the application of different relations. To 
automatize this procedure, we have developed a graph based approach to the SME reduc- 
tion that allows for an efficient implementation and performs a brute-force search for a 
small SME basis. Products of Dirac chains are translated to directed graphs where the 
various structures, like gamma matrices, projection operators, and spinors are represented 
by nodes, and directed edges describe contractions of Lorentz indices and the ordering of 
structures. Algebraic relations then translate to operations on graphs, for instance shrink- 
ing of edges, exchanging or adding of nodes, and result in general in disconnected graphs. 
Since this method can be expressed very efficiently within the framework of graph the- 
ory without the need for computationally expensive algebraic manipulations of lengthy 
expressions, our implementation is capable of testing a huge number of combinations of 
transformations. Typically, the original set of SME is reduced to a basis of several hun- 
dred elements this way. Our variant of the SME reduction will be discussed in more detail 
elsewhere [55]. 

' (1) " (0)* 

As a last step, the products A4n of the NLO SME with the leading-order color 

amplitudes are evaluated, Dirac chains properly contracted, and summations over spins 
and polarizations as well as traces are performed. We translate the resulting expressions 
into C code for an efficient numerical evaluation. 

Altogether, the structure of the final evaluation routines is as follows: each NLO 
diagram/color amplitude contracted with the tree level diagrams is expressed through linear 
combinations of tensor integral coefficients, products of SME, and kinematic invariants. 
The evaluation of tensor coefficients is done as presented in sec. 3, while products of SME 
are computed once per phase-space point and reused. Finally, the sum over squared color 
amplitudes is evaluated. 

It is important to note that, after specifying the desired process and kinematics, no 
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user interaction is required from the point of diagram generation to the final numerical code 
for the cross section evaluation at single phase-space points. All algebraic manipulations 
are performed using FORM while other components like SME and tensor reductions are 
developed in C++. Transparent interfaces, using Python, process input and output between 
the different stages and allow for extensive intermediate checks. The final cross section 
evaluation is made accessible through an automatically generated and flexible C++ interface 
that allows, for instance, the evaluation of single diagrams or color amplitudes interfered 
with the LO contributions, the extraction of divergences, different reduction methods or a 
direct connection with a phase-space generator. 

3 Reduction of tensor integrals 

The integration over the loop momentum in NLO one-loop amplitudes involves N-point 
tensor integrals ^'^^^.../^p that are commonly decomposed into linear combinations of ten- 
sor structures (products of external momenta and metric tensors) with Lorentz invariant 
so-called tensor coefficients T^y The indices {j} encode the rank and composition of the 
corresponding tensor structure. The general strategy for the evaluation of T|J| is their 
reduction to master integrals, usually scalar A^-point functions . In the case of A-point 
tensors with A < 4 the well-known Passarino-Veltman (PV) algorithm [56] can be used, 
while A > 4 coefficients are reduced to linear combinations of four-point tensor integrals. 
Due to numerical instabilities in the vicinity of phase-space points where Gram determi- 
nants become small, alternative reduction techniques exist to produce reliable results. 

Applied methods Our tensor reduction approach combines different methods, allows 
for cross checks between them and ensures numerical stability in an automatized way. For 
A-point functions with A < 4 the following reduction schemes are used: 

• PV reduction [56], 

• reduction with modified Cayley determinants as introduced by Denner and Dittmaier 
(DD) in [57], and 

• expansions around small quantities, like Gram/Cayley determinants and kinematic 
invariants (DDx) developed by the same authors of Ref. [57]. 

In addition, our software is capable of producing multiple precision (MP) reductions with 
help of the qd library [58] , that turn out to be numerically stable already in the framework of 
the PV reduction (MP PV) . Tensor coefficients with five and six external legs are evaluated 
following an approach by Diakonidis et al. [59, 60] that is free of inverse Gram determinants 
(GDF) and therefore numerically stable. 

Our implementation is inspired by [57] and performs the reduction numerically. How- 
ever, the original recursive algorithm is unrolled into an iterative procedure by arranging 
the tensor coefficients in a tree-like structure, which provides fine-grained control over dif- 
ferent aspects of the reduction. We choose the PV reduction for A < 4 in the absence 
of numerical instabilities and the GDF reduction for A > 4 as our standard methods. 
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Based on these reductions, the evaluation tree is constructed for the required set of tensor 
coefficients T^Ji in such a way that the minimal number of evaluations is guaranteed. For 
optimal reuse of intermediate results, coefficients with different mass distributions on inter- 
nal propagators are brought to a standard form with respect to the external momenta and 
internal masses and are treated together. As an example consider a completely massless 
3-point function coefficient that is required for some subdiagrams in the amplitude. The 
same coefficient may appear as a dependency in the reduction of several 4-point functions, 
which can be either massless or contain one massive propagator. The nodes of the evalu- 
ation tree are assigned the default reduction strategy and after initial creation the tree is 
reused for the majority of phase-space points. As already mentioned, this strategy works 
well in large regions of phase space, but becomes numerically unstable if small Gram de- 
terminants in < 4 point coefficients are encountered. In this case, the evaluation tree 
is extended with subtrees for the unstable tensor functions and their dependencies only. 
These newly created subtrees are computed with one of the alternative methods, either DD 
or DDx as needed, or with MP PV, to ensure numerical stability, while all other nodes are 
reduced with the default procedures. While both approaches provide numerically stable 
results, we use the former techniques mainly for cross checks in critical phase-space regions 
while we employ the latter in computations of squared amplitudes. As subtrees are added, 
the reduction program keeps track of the conditions that lead to inconsistencies such that 
the newly created evaluation paths can be reused in future evaluations. 

Numerical stability Detecting numerical instabilities at a single phase-space point is 
in general a non-trivial task without examining the surrounding phase-space domain or 
additional external information. Performing the tensor reduction in d = 4 — 2 e dimensions 
and regulating both ultraviolet and infrared divergences dimensionally, however, offers a 
direct handle on the achieved accuracy. Firstly, the scalar one-loop integrals in terms of 
which the tensor-integral coefficients are reduced have to be known retaining the full pole 
structure. For this task, we use a custom implementation based on QCDLoop [61] for the IR 
poles together with a modified version of LoopTools [62] that allows for multiple precision 
evaluations. During the reduction which is performed on the divergent and finite parts 
separately, UV/IR poles are affected by the same numerical instabilities as the finite part. 
Provided the divergences can be computed for a given tensor coefficient independently in 
a reliable way, a direct comparison can be used to detect a loss of precision^. UV poles 
of tensor coefficients in the minimal subtraction scheme are mass independent and can 
either be obtained with a single analytic reduction or can be taken from the appendix of 
[57]. In contrast, IR divergences occur in certain limits of vanishing kinematic invariants 
and have to be studied in all these cases individually. To this end, we have reduced 
3- and 4-point tensor coefficients up to rank 4 analytically in the various IR divergent 
limits. Using these results a numerically stable library for their computation was carefully 
established and thoroughly cross-checked with the full analytical results for a wide range 
of input parameters. In our reduction, the UV and IR poles of each tensor coefficient are 

^At the cross section level, the same approach to identify numerical instabilities has been successfully 
used in [28, 47, 63]. 
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then checked against this hbrary for every phase-pace point and, in case of inconsistencies, 
alternative reduction methods are employed as described above. Of course, this procedure 
comes with a computational cost, that is mainly due to the multiple reductions for pole and 
finite parts and the evaluation of scalar integrals, while the contribution from the IR pole 
evaluation routines are negligible. However, thanks to the efficient design of the reduction 
algorithm and extensive caching, the run times are competitive with what is reported in 
the literature: computation times in the numerically stable case for all tensor integrals 
required in a mixed massive and massless 2 — > 4 process average at around 20 ms per 
phase-space point on an Intel 17 950 CPU at 3.07GHz. 

4 Characteristics of the automatized approach 

In its current state, our automatized approach is tested for QCD NLO one-loop corrections 
to amplitudes with up to six external particles, of which at least one and at most three 
are (weak or strong) gauge bosons. We have cross checked several 2 — > 3 and 2 — > 4 
processes with the literature. For instance we have reproduced full results for Wbb/ Zbb 
hadroproduction [10-12] and we have checked individual parts of the NLO calculation of 
hard-photon production with heavy quarks (QQ'j for Q = b,t) against an independent 
internal calculation [64]. The 2^4 process ud Wddg was checked at one benchmark 
point against a result obtained in [28]. 

Apart from instabilities in the reduction of tensor coefficients, cancellations in inter- 
mediate expressions of the unrenormalized squared amplitude F in eqn. (2.3) may also 
induce a loss of accuracy in some phase-space regions. In this case, we extend the nu- 
merical precision for both the complete tensor reduction as well as the evaluation of the 
whole contribution to F. Again, the loss of precision is detected by comparison with the 
expected UV divergence structure of F, which is independent of the used renormalization 
scheme and can be obtained with high precision from counterterm contributions. This step 
is computationally most expensive, as a huge number of operations has to be performed 
in slow multiple precision mode both in the tensor reduction and in the evaluation of F. 
Fortunately, the proportion of this type of evaluations is in general relatively small. Com- 
pared to the naive approach where no analysis of instabilities is performed on the tensor 
reduction level, the necessary number of this kind of evaluations is substantially reduced. 

Table 1 gives an overview of the obtained efficiency for the evaluation at 5 • 10^ random 
phase-space points with reasonable cuts, requesting a maximal relative error of 10^^. As 
expected, the evaluation time scales with the number of external particles. Moreover, due 
to a larger basis of SME, amplitudes containing weak couplings compared to for example 
ti'y production are computationally more expensive. An interesting observation, however, 
is the fact that the number of switches to multiple precision evaluations, both within 
the reduction and at the amplitude squared level, do not vary much between processes 
of comparable complexity. Although evaluations in quadruple precision take significantly 
more time with increasing number of external states, the overall evaluation time is governed 
by the numerically stable bulk of phase space. 
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Process 



r, 



s 



Tdq tm/ms tg/ms tq/ms tdq/ms t^"^/ms 



ud ddgW 
ug — )■ bbdW 



qq — >■ jtt 
99 itt ^ 
qq' Whh 
qq Zbb 
gg Zbb 



99.6% 
98.9% 
99.7% 
99.8% 
98.3% 
95.4% 
93.1% 



0.4% 
1.1% 
0.3% 
0.1% 
1.6% 
3.6% 
5.6% 



10.9 10.4 167 1264 

0.1% 17.7 14.4 217 3161 2290 

0.1% 22.5 15.7 233 3314 2706 

1.0% 90.3 37.5 306 4358 5503 

1.3% 95.4 29.7 311 3870 5192 



9.5 8.9 153 1069 



12.0 10.1 182 1972 



Table 1. Benchmarks of the numerically stabilized method applied to various NLO amplitudes for 
the evaluation of 5 • 10^ phase-space points. Ts, and rdq give the ratios of phase-space points 
that required either only standard (double) or also some additional quadruple/double-quadruple 
precision evaluations at the reduction or amplitude-squared level for reliable numerical results. 
gives the mean evaluation time per phase-space point while t^, tq and tdq denote separate mean 
timings for the respective numerical precision. Finally, the mean computation time of both the 
amplitude and tensor reduction in full quadruple precision is given in t^"". The above numbers 
were obtained on an Intel i7 950 CPU at 3.07GHz. 

It is instructive to study the effect of both our approach and different choices of re- 
duction algorithms on the obtained accuracy in the final result. First, the unrenormalized 
squared NLO amplitude of dg — ?• Wbbu is sampled in multiple precision for 5 • 10^ different 
phase-space points and the UV divergent parts are verified to cancel with counterterm 
contributions. The so-obtained results serve as high-precision reference points rref,j. Sub- 
sequently, we compute the squared amplitude for the same set of points using two different 
strategies: 

1. standard reduction of 5- and 6-point tensor coefficients, 

2. GDF reduction of 5- and 6-point tensor coefficients with switches to multiple preci- 
sion both at the n-point tensor integral reduction and amplitude squared level when 
required, requesting a maximal relative error of 10~^ 

and N < 4-point functions are treated with the PV reduction algorithm. Fig. 3 shows the 
distribution of the logarithmic error, 



of the squared amplitude F relative to the reference points Fref for the two different cases. 
Without special handling of numerical instabilities, strategy 1, as expected shows a wide 
distribution of the obtained relative error. Our approach, strategy 2, gives a considerably 
better relative error distribution by reevaluating low precision points. 

5 Result for dg — Wbbu 

For future reference, we provide our new result for the unrenormalized squared amplitude 
of dg Wbbu at NLO at a single phase-space point. The result is normalized to the LO 




(4.1) 
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10000 





Figure 3. Comparison of the obtained numerical precision for the finite part of the process dg — > 
Wbbu at 5 X 10^ phase-space points using standard double precision PV tensor reductions and our 
numerically stable approach. 

cross section in the following way 

. _ r(l - 2e) T 

Svra, r(l + e)r2(l-e) |a^{0)|2' ^ " > 

such that the final result is independent of the strong and weak couplings as well as CKM 
matrix elements. Furthermore, we use 

mw = 80.41 GeV and (5.2) 
mb = 4.62GeV. (5.3) 

for the weak-boson and bottom-quark masses and set all external particles on-shell. 

For dg — )> Wbbu our result with n/ = 4 light and n/^ = 1 heavy-quark flavors at the 
phase-space point of tab. 2 with renormalization scale /x^ = (p^ +Pg)^ reads 

t{dg Wbbu) = -5.6666667 e'^ + 39.342424 e"^ + 292.92493 (5.4) 



6 Conclusions 

We have developed a new automatized approach to the evaluation of one-loop amplitudes in 
terms of Feynman diagrams and applied it to the calculation of the 0{as) virtual corrections 
to qg —7- Wbbq' (and qq' Wbbg). These corrections enter the NNLO calculation of Wbb 
hadroproduction as well as the NLO calculation of both Wbb + j and Wb + j production 
in a fully consistent four-flavor-number scheme. A thorough study of the impact of these 
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H/ 
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P 


P 


Pd 


100.000000000000 








100.000000000000 


Pa 


100.000000000000 








-100.000000000000 


Pu 


14.4169546267975 


-3.59819144566031 


6.52544251406004 


-12.3418069595668 


Pb 


53.6542637065835 


-16.9076522158373 


-49.1575349754512 


12.4540622120327 


Pb 


25.2318438952597 


-17.2383739318242 


-15.9080092164594 


8.06692341047065 


Pw 


106.696937771359 


37.7442175933219 


58.5401016778506 


-8.17917866293656 



Table 2. Phase-space point used for dg Wbbu 

corrections in both previous cases as well as the application of the method developed in 
this paper to other processes will be the subject of future publications. 
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